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ABSTRACT 

We present a catalogue of 1602 white dwarf-main sequence binaries (WDMS) from the 
spectroscopic Sloan Digital Sky Survey Data Release 6 (SDSS DR6). Among these we 
identify 440 as new WDMS binaries. We select WDMS binary candidates from tem- 
plate fitting all 1.27 million DR6 spectra, using combined constraints in both % 2 and 
signal-to-noise ratio. In addition, we use Galaxy Evolution Explorer (GALEX) and 
UKIRT Infrared Sky Survey (UKIDSS) magnitudes to search for objects in which one 
of the two components dominates the SDSS spectrum. We use a decomposition/fitting 
technique to measure the effective temperatures, surface gravities, masses and dis- 
tances to the white dwarfs, as well as the spectral types and distances to the com- 
panions in our catalogue. Distributions and density maps obtained from these stellar 
parameters are then used to study both the general properties and the selection ef- 
fects of WDMS binaries in SDSS. A comparison between the distances measured to 
the white dwarfs and the main sequence companions shows d scc > c? w d for ~l/5 of the 
systems, a tendency found already in our previous work. The hypothesis that mag- 
netic activity raises the temperature of the inter-spot regions in active stars that are 
heavily covered by cool spots, leading to a bluer optical colour compared to inactive 
stars, remains the best explanation for this behaviour. We also make use of SDSS- 
GALEX-UKIDSS magnitudes to investigate the distribution of WDMS binaries, as 
well as their white dwarf effective temperatures and companion star spectral types, 
in ultraviolet to infrared colour space. We show that WDMS binaries can be very 
efficiently separated from single main sequence stars and white dwarfs when using a 
combined ultraviolet, optical, and infrared colour selection. Finally, we also provide 
radial velocities for 1068 systems measured from the Nai AA 8183.27,8194.81 absorp- 
tion doublet and/or the Ha emission line. Among the systems with multiple SDSS 
spectroscopy, we find five new systems exhibiting significant radial velocity variations, 
identifying them as post-common-envelope binary candidates. 

Key words: Binaries: spectroscopic - stars:low-mass - stars: white dwarfs - binaries: 
close - stars: post-AGB - stars: evolution variables 



1 INTRODUCTION 

Binaries containing a white dwarf primary plus a main 
sequence companion were initially main sequence binaries 
in which the more massive star evolved through the gi- 
ant phase and became a white dwarf. In the majority of 
cases the initial separation of the main sequence binary 
is wide enough to allow the evolution of both stars as if 



they were single. A small fraction is believed to undergo a 
phase of dynamically unstable mass transfer once the more 
massive sta r is on the giant branch or the asymptotic gi- 
ant branch l|Webbinklll98i Ide Koo]|[l992l IWillems fc Kolbl 
2004). As a consequence of this mass-transfer the envelope 
of the giant will engulf its core and the companion star, i.e. 
the system is entering a common envelope phase (CE, e.g. 
iLivio fc Sokeall988l : llben fc Liviolll995l : iTaam fc Sandquistl 
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l200d ; IWebbinkllioOsh . Friction inside this envelope causes a 
rapid decrease of the binary separation. Henceforth orbital 
energy and angular momentum are extracted from the bi- 
nary orbit and lead to the ejection of the envelope, exposing 
a post-common-envelope binary (PCEB). As a consequence 
the orbital period distribution of WDMS binaries is clearly 
bi-modal, with PCEBs concentrated at short orbital peri- 
ods, and wide WDMS binaries ( non-PCEBs) at long orbital 
periods (|Willems fc Kolbl 12004! ) . After the ejection of the 
envelope, close WDMS binaries evolve to shorter orbital pe- 
riods through angular momentum loss via magnetic braking 
and/or through the emission of gravitational waves. WDMS 
binaries include progenitors of a wide range of astronomical 
objects such as cataclysmic variables and super-soft X-ray 
sources, with some of those objects likely evolving at later 
stages into type la supernova. 

Population synthesis models have been developed 
for a variety of bin a ry st ar s undergoing CE evo- 
lution (iDewi fc Taurisl 
iPolitano fc Weilerl 2006 , 



20001. iNelemans fc Toutl 12005 
20071 ; iDavis et al.l |200S| . |2009[ ) 



However, the theoretical understanding of both CE evolu- 
tion and magnet ic braking is currently poor ly constrained 
by observations IjSchreiber fc Gansickdl2003l ). and progress 
on this front is most likely to arise from the analysis of a 
large sample of PCEBs that are well-understood in terms 
of their stellar components. WDMS binaries appear most 
promising in that respect, as their st ellar components 
are relatively s i mple, and the SPSS ( York et all |2000| ; 
IStoughton et all |2002| ; lAbazaiian et alj l2009h offers the 
possibility to dramatically increase the number of WDMS 
binaries available for detailed follow-up studies. Up to date 
there e xist four compil a tions of SDSS WDMS binaries, 
namely lEisenstein et al.l (120061) (obtained a s part of their 
search of w hite dwarfs) , Silvestri et al.l J2007h ( which 
contains also iRavmond et al.1 (|2003h : ISilvestri et all (120061 ) 
as a s ubset, and was cla i med t o be complete for the SDSS 
DR5), lAugusteiin et al.1 i|200Sl ) (obtained from SDSS DR5 
using colour cuts plus proper motions), and iHeller et al.l 
(2009) (obtained from a search for sdB stars in the SDSS 
DR6). 

Here we initiate a comprehensive study to compile a 
master sample of spectroscopic WDMS binaries from SDSS. 
In this first pu blication we make use of the S DSS spec- 
troscopic DR6 (lAdelman-McCarthv et ail [2008) to create 
a catalogue of 1602 WDMS binaries and candidates that 
were serendipitously observed. This list is not only a su- 
perset of all those previous SDSS WDMS binary compila- 
tions, but also includes 440 new WDMS binaries. In ad- 
dition we provide a coherent analysis of the system pa- 
rameters of both stellar components, as well as an exten- 
sion to GALEX/UKIDSS colours to study the properties 
of SDSS WDMS binaries. In a parallel effort, and as part 
of SEGUE (the SDSS Extension for Galactic Understand- 
ing and Exploration), some of us carried out a dedicated 
program to id entify ~ 300 WD MS binaries containing cold 
white dwarfs (|Schreiber et al.ll2007l ). a population clearly 
underrepresented in previous samples of WDMS binaries. A 
detailed analysis of the SEGUE population of WDMS is in 
preparation (Schreiber et al. 2009). Finally, we plan to sum- 
marise our efforts by presenting the WDMS binary conten t 
of the final SDSS data release (DR7. lAbazaiian et all l2009). 



This last paper will be accompanied by a public online data 
base of WDMS binaries. 

The final master sample of "all" spectroscopic SDSS 
WDMS binaries will form a superb database for fu- 
ture follow-up studies of WDMS binaries. In partic- 
ular, analysing the fraction and the characteristics of 
the PCEBs among the WDMS binaries may provide 
strong con straints on current theor i es of compact binary 
evolu t ion dRebassa-Mansergas et al.l 12007 \_ Schreiber et al.1 



20081: iRebassa-Mansergas et al.l |200S|; iPvrzas et al.l 120091 ; 
Nebot Gomez- Moran et alj|2009l ; ISchwope et al.ll2009l ). 

The structure of the paper is as follows. In Sect. [2] we 
present our method of identifying WDMS binaries, and es- 
timate the completeness of the sample in Sect. [3] In Sect.f?] 
we provide our final catalogue. Using a spectral decomposi- 
tion/model atmosphere analysis, we derive white dwarf ef- 
fective temperatures, surface gravities, masses, companion 
star spectral types, and distances in Sect .[5] We compare 
these stellar parameters to those obtained in other stud- 
ies in Sect.[7] and study the selection effects of WDMS 
binaries within SDSS in Sect. [8] We provide then colour- 
colour diagrams and colour-colour cuts for WDMS binaries 
in Sect.[5]and Sect. 1101 respectively. We finally measure the 
NaiAA 8183.27,8194.81 absorption doublet and/or the Ha 
emission radial velocities in Sec. llll and summarise our re- 
sults in Sect. 1121 



2 IDENTIFICATION OF WDMS BINARIES IN 
SDSS 

2.1 Computational method 

We have developed a procedure based on x 2 template 
fitting in order to automatically identify WDMS binary 
candidates from the SDSS DR 6 spectroscopic data base 
l|Adelman-McCarthv et ail 120081 ) . Our initial template set 
consisted of se veral dozen SDSS s pectra of con firmed WDMS 
binari es from lEisenstein et al.l (|2006T > and ISilvestri et al.l 
(2007). These spectra were chosen to sample a broad range 
in white dwarf temperatures and subtypes (DA, DB, and 
DC) as well as companion star spectral types, and to be 
of high signal-to-noise ratio (S/N). A set of representa- 
tive templates is shown in Fig.[T] In addition, we compiled 
a set of 17 sing le DA w hite dwarf template spectra from 
lEisenstein et all 's (|2006l ) list, covering the entire observed 
range of T e g and log g, as well as the M0-M9rBochansk i et al.l 
l|2007l ) M-dwarf templates. 

Each of these WDMS binary, white dwarf, and M-dwarf 
templates were then fitted to the full 1.27 million spectra in 
DR6. In this process, the template spectrum was normalised 
to the SDSS spectrum under scrutiny, and a reduced \ 2 was 
calculated using the flux errors of the two spectra added 
in quadrature. In practice, our fitting procedure produced 
for each of the WDMS binary, white dwarf, and M-dwarf 
templates a list of spectrum identifier (MJD-PLT-FIB), S/N 
of the spectrum, and \ 2 f° r all SDSS DR6 spectra. For each 
of the templates, we plotted x 2 as a function of S/N of the 
target spectrum (see Fig[5]), and defined a power law relation 
Xmax — a x S/N 6 . We considered any spectrum with 
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Figure 1. Six examples of previously known WDMS binaries used in this work as WDMS binary templates. SDSS names and MJD- 
PLT-FIB identifiers are indicated for each of them. White dwarf effective temperatures and spectral type of the companions (see Sec.[5) 
are also indicated in each panel. 



as a candidate WDMS binary, white dwarf, or M-dwarf (de- 
pending on the current template). The S/N-dependent form 
of Xmax accounts for the increase of x 2 f° r higher values of 
S/N. This constraint had to be defined individually for each 
of the templates, as the different spectral shapes resulted 
in a large spread of % 2 distributions. The (x 2 , S/N) planes 
obtained from fitting the SDSS spectra are shown for two 
different WDMS binary templates in Fig.[U 

After a first run through the DR6 spectra, we comple- 
mented the template set with the spectra of a number of 
newly identified WDMS binaries, and re-run the fitting for 
those new templates again. That process was repeated un- 
til no new WDMS binary candidates were found - at which 
point we had used a total of 163 different WDMS binary 
template spectra. 

Even though the above method efficiently identifies 
WDMS binary candidates among the spectra in DR6, the 
choice of Xmax alone does not avoid completely the filter- 
ing of other astronomical objects, such as quasars, main se- 
quence stars, and galaxies. In addition, for templates that 
are dominated by the white dwarf (M-dwarf), the list of 
candidates will unavoidably contain a substantial number 



of single white dwarfs (M-dwarfs). Hence we first visually 
inspected all WDMS binary candidates, as well as the white 
dwarf and M-type star subsamples (a total of ~ 70000 spec- 
tra), and removed those objects that were not of our in- 
terest, i.e. neither WDMS binary, white dwarf, or M-dwarf 
candidates. The final result of the template fitting were a 
list of 1491 WDMS binary candidates, 8368 single white 
dwarf candidates, and 15379 single M-dwarf candidates. It 
is worth mentioning that we found here 36 WDMS binaries 
that were observed only/first by SEGUE, and consequently 
decided to include them in the SEGUE list of WDMS bina- 
ries (Schreiber et al. 2009, in preparation). 

2.2 Red and blue excess in SDSS spectra: help 
from GALEX and UKIDSS 

While the template fitting proved to be a robust method 
to find WDMS binaries in which both stellar components 
contribute clearly visible amounts of flux, the procedure is 
prone to mis-classify white dwarf dominated WDMS bina- 
ries as single white dwarfs, and M-dwarf dominated WDMS 
binaries as single M-dwarfs. We therefore decided to probe 
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more specifically for the presence of excess flux at the red 
(blue) end of the SDSS spectra in objects classified initially 
as single white dwarfs ( M- dwarfs ) . 

For the search of red flux excess in single white dwarf 
candidates, we fitted synthetic w hite dwarf spe c tra c om- 
puted with the code described by iKoester et al.l l|2005h to 
the SDSS spectra, and then calculated the reduced \ 2 over 
the wavelength ranges 4000-7000 A (x?) and 7000-9000 A 
(Xi)- Objects with Xi/Xb > 1.5 were "promoted" from single 
white dwarf candidates to WDMS binary candidates. 

The search for blue flux excess proceeded in an analo- 
gous fashion for the single M-dwarf candidates, only that 
we used the set of hig h S/N M-dwarf templates from 
iRebassa-Mansergas et aL (|2007f ) instead of model spectra, 
and calculated the reduced \ 2 over the wavelength ranges 
4000-5000 A and 7000-9000 A. Objects with xt/Xr > 1-5 
were "promoted" from single M-dwarf candidates to WDMS 
binary candidates. 

On the left and right top panels of Fig. [3] we 
show the SDSS spectra (black line) and SDSS mag- 
nitudes (red dots) of SDSSJ 132925.21+123025.5 and 
SDSSJ 131928.80+580634.2, along with the best-fit white 
dwarf model and M-dwarf template (red lines, middle pan- 
els). These two objects were initially classified by our tem- 
plate fitting procedure as single white dwarf and single M 
star candidates, respectively, but "promoted" to WDMS bi- 
nary candidates by the flux excess measurement as described 
above. The flux excess is more obvious when plotting F v 
(middle left panel) instead of Fx (top left panel). However, 
in several cases, the detection of blue or red flux excess is 
rather marginal. 

As a final step in our search for WDMS binaries, we 
have cross-correlated o ur entire list of WDMS binary can- 
didat es with GALEX |Martin et al.l 120051 ; iMorrissev et all 
l2005h DR4, providing near- and far-ultraviolet ( nuv, fuv) 
magnitudes, and with the PR 4 of UKIDSS (|Dve et al l 
l2006l ; iHewett et al-llioOrj ; lLawrence et al.l l2007h . providing 
infrared yJHK magnitudes. We then inspected the ob- 
served ultraviolet-optical spectral energy distribution of all 
secondary star dominated WDMS binary candidates, and 
the optical-infrared spectral energy distribution of all white- 
dwarf dominated WDMS binary candidates. Objects where 
a clear ultraviolet or infrared excess was detected were then 
included in our WDMS binary sample. 

For SDSSJ 132925.21+123025.5, the UKIDSS magni- 
tudes unambiguously confirm the existence of a low-mass 
companion (bottom left panel of Fig.[3J). Similarly, the ultra- 
violet GALEX magnitudes clearly confirm the presence of a 
white dwarf primary in SDSSJ 131928.80+580634.2 (bottom 
right panel of Fig. [3} . 

2.3 SDSS images 

As a final check on the nature of the WDMS binary candi- 
dates, we inspected their SDSS DR6 images for morpholog- 
ical problems, and found primarily two types of issues. 

Firstly, single white dwarfs (M-dwarfs) may occasion- 
ally be located close to very bright M-dwarfs (white dwarfs 
or A-stars) causing scattered light to enter the spectroscopic 
fibre resulting in an (apparent) two-component spectrum. 
A spectacular example is SDSSJ 073531.86+315015.2 (left 
panels of Fig.[4j : the SDSS spectrum clearly exhibits an M- 



dwarf at red wavelengths and a blue component with strong 
Balmer lines in the blue - however, the SDSS image reveals 
that this is a single M-dwarf at a distance of 12 arc-min of 
Castor A/B - two V = 1.88 and V = 2.96 A-stars. The 
SDSS magnitudes (red dots) are superimposed to the SDSS 
spectrum (black) and are consistent with a single red star. 
Single M-dwarfs are also likely to be found superimposed 
with single early-type stars in the same figure. An exam- 
ple is SDSSJ 005827.24+005642.6 (see middle top panel of 
Fig-El. At first glance one could be tempted to consider it 
a spatially resolved WDMS binary pair. However, the SDSS 
spectrum (middle bottom panel of the same figure) shows 
the typical Balmer lines of an early F star in the blue, while 
at redder wavelengths the typical spectral features of a low- 
mass star can be seen. The large difference in absolute mag- 
nitudes between the two spectral types implies that these 
two stars are a chance superposition along the line of sight, 
and not a physical binary. 

Secondly, SDSS images can help identifying WDMS 
binaries among our sample that are spatially resolved, 
but close enough that flux from both stars will enter 
into the spectroscopic fibre. In such cases, the SDSS 
magnitudes are often discrepant with the flux-calibrated 
SDSS spectrum, and/or have large errors as consequence 
from the deblending applied by the photometric pipeline. 
Figure. [4] (top right panel) shows the SDSS image of 
SDSSJ 025306.37+001329.2, which reveals a spatially re- 
solved pair of red and blue stars. The SDSS spectrum of 
SDSSJ 025306.37+001329.2 contains the typical signatures 
of a WDMS binary, i.e. broad Balmer lines from the white 
dwarf and TiO absorption bands from the M-dwarf, how- 
ever, the errors on the SDSS magnitudes are untypically 
large, and do not match well the flux calibrated SDSS spec- 
trum. 

2.4 Cross-checks with previous WDMS binary 
catalogues 

A total number of 1491 WDMS binary candidates were iden- 
tified in Sect. l2.ll From the analysis carried out in Sec. 12.21 
and Sec.[23] we have identified 94 and 89 WDMS bina- 
ries by their blue and red excess respectively, while 115 
WDMS binary candidates were removed after inspecting 
their SDSS images. This increased the number of systems 
to 1559 WDMS candidates in the spectroscopic SDSS DR6 
data base. 

In order to evaluate the efficiency of our proce- 
dure we compared our results to five previously pub- 
lished lists of WDM S bina rie s from SPSS, namel y 
van den Besselaar et alj (|2005h . lEisenstein et al.l lj2006h . 
Silvestri et al.l (120071) (whi ch includes the sys tems from 
Raymond et""ai1 (120031) and ISilvestri et all d2006f> as a sub- 
set). I Augusteiin et all (j2008i ). and lHeller et all (I2009T ). This 
comparison is summarised in TableQ] 

The 15 WPMS binaries containing PB (13) and DC (2) 
white dwarfs presented by Ivan den Besselaar et alj (|2005l ) 
have been successfully identified as WDMS binaries with 
our template fitting algorithm. 

The WDMS binary sample presented by 
lEisenstein et al.l (|2006h is almost entirely contained in 
the W DMS binary catalogue provided by ISilvestri et al.l 
l|2007l ) (see below for a detailed comparison). Only ten 
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Figure 2. x 2 -S/N plane obtained fitting two of our WDMS binary templates (SDSS J103121.97+202315.1, left and SDSS J204431.44- 
061440.2, right) to the entire SDSS spectra database. Objects falling in the area defined by Xspcc < 1 X S/N 6 were considered WDMS 
binary candidates. Left panel: Xspcc < 0.001 X S/N 31 . Right panel: Xspcc < 0.0001 X S/N 3 ' 05 . WDMS binary candidates are shown in 
magenta, in red the equation Xspcc < a x S/N 6 denned for each template. 
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Figure 3. Top left: SDSS spectrum of SDSSJ 132925.21+123025.5, a WDMS binary candidate initially catalogued as single DA 
white dwarf. The red dots represent the SDSS magnitudes. Middle left panel: the best white dwarf model fit is superimposed in 
red, unambiguously identifying the red excess of the binary. Bottom left panel: SDSS and UKIDSS magnitudes superimposed to the 
SDSS spectrum. Again, the UKIDSS magnitudes clearly show the presence of a low-mass companion. Top right panel: the same for 
SDSSJ 131928.80+580634.2, an initially catalogued early M-type star. Middle and bottom right panels: the best M-type fit and the 
ultraviolet GALEX magnitudes clearly confirm the presence of a white dwarf primary, respectively. The red and blue straight lines 
represent the white dwarf solutions (red for the hot, blue for the cold) obtained from decomposing/fitting the spectrum (see Sect[5jl. 
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Table 1. Comparison of the WDMS binaries identified here with 
those provided in previous publications. Nwdms gives the num- 
ber of objects listed, N;(j en t the number of WDMS binaries we 
have identified within the list, N remov the number of objects we 
do not consider being WDMS binaries, and Ni ost the number of 
systems that we have missed (the percentage is given in brackets) . 



Publication 


Nwdms 


Nident 


Nremov 


Ni os t 


van den Bessclaar ct al. (2005) 


15 


15 




- (0%) 


Silvestri et al. ("20071 


1225 


996 


204 


25 (2%) 


Augusteiin et al. (2008) 


130 


110 


16 


4 (3%) 


Heller et al. (2009) 


636 


554 


81 


1 (0.1%) 



Tabl e 2. Updated classification of the 204 objects from 
ISilvestri et"aTTj2007f) which are not considered as WDMS binaries 
by us. The complete table can be found in the electronic edition 
of the paper. 



Object 



Classification 



SDSSJ001324.33-085021.4 
SDSSJ003839.36+260258.5 
SDSSJ005714.52-000755.8 
SDSSJ005827.24+005642.6 
SDSSJ012516. 98-010944.2 



M star + ? 
G star 
M star 

MS+MS superposition 
M star 



We identify nine quasars, three cataclysmic variables, one F 
star, one M-dwarf, one DA white dwarf, and one DC white 
dwarf. The four remaining WDMS binaries have been over- 
looked by our method because they contain rather cold white 
dwarfs and the spectra are dominated by the companions. 
We added th e four systems to o ur sample. 

: Recently iHeller et al.l (|2009l ) discussed the properties of 
636 colour selected WDMS binaries from SDSS DR6. Most 

"of their systems had been identified previously. Comparing 
their sample with our catalogue we find 82 objects miss- 
ing in our list. Inspecting the SDSS images and spectra (if 
available) of these 82 systems we find 81 of them not being 

-WDMS binaries. Two examples are the low-redshift galax y 
SDSSJ 122953.46+473150.3 (U7639, iMakarova et"allll998h . 
and again the z = 0.21 quasar SDSS J 032428.78-004613.8, 
also included in lSilvestri et al] (|2007l ). An updated classifi- 
cation for these 81 objects is provided in Table[3] whilst the 
complete table is provided in the e lectronic edition o f this 
paper. Only one system analysed bv lHeller et al. (2009]) but 
not in our sample turned out to be a WDMS binary, and we 
added it to our catalogue. 

In summary, comparing the sample of WDMS binaries 
obtained with our template fitting procedure with previ- 
ously published catalogues of WDMS binaries shows that 
our method represents a robust and efficient tool to identify 
WDMS binaries (Tabled). 



objects of lEisenstein et al. | d2006h with avail able SDSS 
spectra are not listed in ISilvestri et al.l (2007), and our 
algorithm has successfully identified these ten systems as 
WD MS binaries. 

ISilvestri et all (|2007l ) claim that their catalogue con- 
tains 1253 objects but in fact only 1228 spectra (correspond- 
ing to 1225 objects) are listed in the electronic edition of 
their paper, and 996 systems of those appear in our cat- 
alogue. One is a SEGUE object and will be included in 
Schreiber et al. (2009, in preparation). The vast majority 
(i.e. 204) of the remaining 229 objects are not classified as 
WDMS binaries by us, either because of their spectroscopic 
appearance, or because of morphological problems in the 
SDSS images. Two examples of mis-classified WDMS bi- 
naries are the z = 0.21 quasar SDSS J032428. 78-004613.8 
and the z = 0.11 galaxy SDSS J114334.70+455134.2. An 
updated classification of the 204 objects is given in Tabled 
whilst the full table is available in the electronic edition of 
this paper. T he spectra of the 25 genuine WDMS binaries 
identified by ISilvestri et al.l (2007) but overlooked by our 
template fitting algorithm, are dominated by the emission 
of one of the stellar components. We added these 25 WDMS 
bin aries to our sampl e . 

lAugusteiin et al.l l|2008l) developed a WDMS binary 
identification algorithm based on SDS S imaging and the 
prope r motion catalogue provided by I Gould fc Kollmeierl 
lj2004) . This way they identified 651 WDMS binary can- 
didates with SDSS ugriz photometry, of which 95 are con- 
tained in the SDSS DR5 spectroscopic data base. Cross- 
correlating their full list of 651 objects against the DR6 
spectroscopic data base we find 176 spectra for 130 objects. 
All but 20 of these objects were in our list. Again, the ma- 
jority of those systems (16) that have not been identified 
by our template fitting procedure are not WDMS binaries. 



3 COMPLETENESS OF THE SAMPLE 

In the previous section we showed that our selection method 
successfully recovered the vast majority of previously iden- 
tified SDSS WDMS binaries. In this section, we investi- 
gate both the internal completeness of our catalogue, i.e. 
the fraction of WDMS binaries contained in the DR6 spec- 
troscopic data base that has been identified by our algo- 
rithm, and the external completeness, i.e. the fraction of 
point sources for which SDSS ugriz photometry is available 
that has been spectroscopically observed by SDSS. A de- 
tailed study of the external (spectroscopic) completeness of 
WDMS binaries within SDSS is currently underway, and we 
use this paper to provide preliminary results on this issue. 
To that end, we first define three regions in SDSS colour- 
colour planes: t he WDMS bina r y exc lusion region that has 
been defined by iRichards et al.l ( 20021 ) as part of the quasar 
fibre allocation algorithm, and two small rectangular boxes 
in the (g — r,r — i) plane that sample a significant popula- 
tion of WDMS binaries in colour space. Particularly, these 
two boxes were defined to test how the external complete- 
ness varies between inside and outside the WDMS exclusion 
box. The (u~ g,g — r), (g — r,r — i), and (r — i, i — z) colour- 
colour diagrams shown in Fig. [5] illustrate the locations of 
these regions (WDMS, box 1, box 2). The WDMS binaries 
in our sample excluding those systems that are classified 
as candidates only (see Sect.|4)| are shown in black, stellar 
sources are in gray, and quasars in light gray. Highlighted in 
light blue are the white dwarf (WD), A-star (A) and WDM S 
binary exclusion regions defined by IRichards et all i|2002l ). 
Both the A-star and the WDMS binary region in the colour 
planes are combined with a logical "and" . Box 1 and box 2 
are defined by (-0.5 < g - r < -0.2, < r - i < 0.2) and 
(—0.3 < g — r < 0.1, 0.3 < r — i < 0.5) respectively. 
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Figure 4. SDSS images of three WDMS binary candidates. Top left panel: SDSS image of SDSSJ 073531.86+315015.2 (30' X 30'), a 
single red star initially considered as a WDMS binary candidate. Bottom left panel: SDSS magnitudes (red dots) and spectrum (black 
line) of the same system. The light from the saturated bright star (Castor A/B) is also dispersed in the spectrum. The magnitudes are 
consistent with a single red star. Top middle panel: SDSS image of SDSSJ 005827.24+005642.6 (1' X 1'). The image suggests a resolved 
WDMS binary pair. Middle bottom panel: the detection of the Balmer lines typical of an early F star in the blue, together with the typical 
spectral features of a low-mass main sequence star in the red (black solid line), indicate that these are two single stars superimposed in 
the same image rather than a resolved WDMS binary pair. SDSS magnitudes are indicated with red dots, and are consistent with those 
of a low-mass star. Top right panel: SDSS image of SDSSJ 025306.37+001329.2 (1' X 1'), a resolved WDMS binary in our sample. Bottom 
right panel: SDSS magnitudes (red dots) and spectrum (black line) of the same system. Whilst the SDSS spectrum clearly shows both 
components, the SDSS magnitude errors are untypically large due to unsuccessful dcblcnding of the close pair. 



To determine the completenesses we visually classified 
all SDSS DR6 spectr a of point-sources w ith g ^ 20, and used 
the casjobs interface dLi fc Thakaj|2008l ) Fl to the SDSS data 
base to determine the number of photometric point sources 
with clean photometry in each of the regions defined above. 
The fraction of WDMS binaries identified by the template 
matching routine among the total number of spectroscop- 
ically identified WDMS binaries inside the selected region 
gives an estimate of the internal completeness. The external 
completeness is simply given by the fraction of spectroscop- 
ically observed point sources inside a given colour-colour re- 
gion. In order to evaluate the impact of the brightness limit 
appli ed in the quasar selection algorithm (Ric hards et al.l 
|2002| ). i.e. a de-reddened i < 19.1, we additionally performed 
the external completeness analysis for the ranges i ^ 19.1 
and 19.1 < i ^ 20.1 outside the WDMS exclusion box, i.e. 



Table 3. Updat ed classification of the 82 objects from 
iHeller et all (120091') which are not considered as WDMS binaries 
by us. The complete table can be found in the electronic edition 
of the paper. 



Object 



Classification 



SDSSJ013007.13+002635.3 
SDSSJ020538.10+005835.3 
SDSSJ023906.69+002916.6 
SDSSJ025622.18+330944.8 
SDSSJ032428.78-004613.8 



carbon star 
no available spectra, WDMS? 

DA 
O star? 
quasar 



for box 1 and box 2. The results are given in Table|4] and 
can be summarised as follows. 



http: / /casjobs. sdss.org/Cas Jobs/ 



• SDSS DR6 cont ains 8002 photometri c point sources 
with g < 20 inside the lRichards et all (|2002h WDMS binary 
exclusion box. SDSS spectra are available for 575 (or 7% 
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Figure 5. From left to right the (u — g,g — r), (g — r,r — i) and (r — i.i — z) colour diagrams for the WDMS binaries in our sample (black; 
excluding WDMS binary candidates), ste l lar so urces (dark ray), and quasars (light gray). The white dwarf, A-star and WDMS binary 
exclusion boxes defined bv lRichards et aL I ll2002l) are shaded in lig ht blue. The regions within the two red boxes in the (g — r,r — i) plane, 
and the WDMS binary exclusion region of iRichards et all [|2002l ) in the (g — r,r — i) and (r — i,i — z) planes are used to estimate the 
completeness of our WDMS binary sample (see Sect. [3) ■ 



of the point sources) systems, and we visually classify 389 
of these as WDMS binaries. The remaining 186 objects are 
mostly single M stars but we also identify four cataclysmic 
variables and eight quasars. The small external spectro- 
scopic completeness of SDSS DR6 (7%) in the WDMS binary 
box as well as the large number of WDMS binaries (67%) 
among the spectroscopically observed objects is not surpris- 
ing, as this region h as been explicitly excl uded in the quasar 
selection algorithm ( Ri chards et al.|[2002l ). All but 12 of the 
389 WDMS binaries have been successfully identified by our 
automated template matching algorithm, which gives an in- 
ternal completeness of 97%. Among the 12 missing objects 
nine contain cool (probably DC) white dwarfs, one contains 
a clear DA white dwarf, one a clear DC white dwarf primary, 
and one we identify as a lo w-accretion-rate pol ar, LARP 
(SDSS J204827.90+005008.9. ISchmidt et al1l2005l ). In all the 
12 cases the spectra are dominated by the emission of the 
companion stars. 

• Box 1 contains 708 point sources with g < 20 in DR6 
and for 247 of these at least one SDSS spectrum is avail- 
able, corresponding to an external completeness of 35%. As 
expected, we find the external completeness to significantly 
change at the i-magnitude limit incorporated in the quasar 
search algorithm. While it increases to ~ 60% for i ^ 19.1 
it drops to ~ 15% for 19.1 < i ^ 20.1. 

Among the 247 objects are 67 (27%) WDMS binaries. The 
remaining systems are mainly quasars, a few single white 
dwarfs, early-type main sequence stars, and two cataclysmic 
variables. All the 67 WDMS binaries in this box have been 
identified by our search algorithm (Sect.[2]| equivalent to an 
internal completeness of 100%. 

• SDSS DR 6 contains 6689 SDSS point sources in box 2 
and for 2280 of them SDSS spectroscopy is available, which 
gives an external completeness of 34% using g < 20. 

As for box 1, the quasar limit on the reddening corrected 
i-magnitude of 19.1 is very much affecting the external com- 
pleteness. It increases significantly to ~ 60% for i ^ 19.1 
but decreases to ~ 14% for objects fainter than the quasar 



Table 4. Wc give here the number of SDSS point sources N p , the 
number of available SDSS spectra N spcc , the number of spectro- 
scopic WDMS binaries N\vdmSj and the number of WDMS bina- 
ries we have identified N;,j eIlt for each of the WDMS, box 1 and 
box 2 regions defined in Sect. [3] The external completeness is esti- 
mated as N S pec/Np, the internal completeness as Njdent /NwDMS ■ 





Reg. 


N p 


Nspcc 


NwDMS 


Nident 


Ext.C. 


Int.C. 


9 <20 


WDMS 


8002 


575 


389 


377 


7% 


97% 




Box 1 


708 


247 


67 


67 


35% 


100% 




Box 2 


6689 


2280 


135 


131 


34% 


97% 


i < 19.1 


Box 1 


313 


188 






60% 






Box 2 


2822 


1672 






60% 




i < 20.1 


Box 1 


389 


59 






15% 




i > 19.1 


Box 2 


13110 


1770 






14% 





limit, i.e. for 19.1 < i ^ 20.1. At first glance, one might con- 
sider the spectroscopic completeness of ~60% for i ^ 19.1 
to be surprisingly low compared to the ~ 95% obtained fo r 
the quasar selection algorithm (IVanden Berk et al.l 120051 ) . 
However, one should keep in mind that "outside the exclu- 
sion boxes" is not equivalent to inside the (rather complex) 
quasar selection. In addition, in DR6 the area of the foot- 
print followed-up spectroscopically was only ~ 82% of the 
imaging area, implying a systematically lower spectroscopic 
completeness than expected for DR7. As mentioned earlier, a 
more complete analysis of the completeness of SDSS WDMS 
binaries is underway and will be presented in a subsequent 
publication. 

From the 2280 spectra that are available in box 2 only 
135 (6%) are WDMS binaries. The vast majority of the re- 
maining objects are quasars but we also identify some few 
single main sequence stars and one cataclysmic variable. All 
but four of the 135 WDMS binaries have been found by 
our template fitting method, corresponding to an internal 
completeness of 97%. The four systems that have not been 
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identified by our algorithm are SDSS J110539. 77+250628. 6, 
which is in fact the semi-detached magnetic cataclysmic 
variable STLMi observed during a deep low state; 
SDSS J124959. 76+035726.6, a typical WDMS binary pre- 
viously listed as a cataclysmic variable candidate by 
ISzkodv et all (|2004h : SDSS J150954.40+243449.3, which has 
a broken SDSS spectrum; and SDSS J204218. 52-065638.4, a 
spatially resolved main sequence dominated WDMS binary. 

Three main conclusions can be drawn from the analy- 
sis carried out in Fig. [5] Firstly, the spectroscopic complete- 
ness in SDSS is much larger in quasar dominated colour- 
colour space, i.e. 34% in box 2, co mpared to only 7% i n 
the WDMS binary exclusion box of iRichards et all (2002). 
The low completeness in the WDMS binary exclusion box, 
combined with the high fraction of WDMS binaries among 
all objects in this region (67%) implies that the number of 
SDSS WDMS binaries could be dramatically increased by 
additional follow-up spectroscopy of point sources located 
in this box. Since the spectroscopic completeness in the 
WDMS binary exclusion box is 7%, and only 389 WDMS 
binaries benefit from spectra, this implies then that SDSS 
did not target ^5500 WDMS binaries for follow-up spec- 
troscopy within this region. As partners of SEGUE some 
of us performed s uch a program and id entified ~ 300 new 
WDMS binaries jSchreiber et al.ll2007h . Secondly, outside 
the exclusion boxes the external completeness drops signif- 
icantly from ~ 60% to ~ 15% at the i-magnitude limit of 
19.1 implemented in the quasar search algorithm. Finally, 
and most importantly, only 16 WDMS binaries (four in box 2 
and 12 in the WDMS binary exclusion box), have not been 
identified by our template fitting routine. Assuming that the 
three analysed colour boxes are representative for the entire 
WDMS binary bridge, we derive an internal completeness of 
our catalogue of ~ 98%. Virtually all the systems our search 
algorithm failed to identify are dominated by the emission 
of the secondary star. Such systems are therefore expected 
to form the remaining ~ 2% of WDMS binaries contained 
in the SDSS DR6 data base. 

The 16 WDMS binaries that we previously overlooked 
have been added to our sample that now forms our final 
SDSS DR6 WDMS catalogue, as described in the following 
section. 



4 THE FINAL CATALOGUE 

From the analysis described in Sect .[2] and Sect .[3] a total 
number of 1602 WDMS binaries have been identified. These 
systems form our final catalogue of SDSS DR6 WDMS bi- 
naries and WDMS binary candidates. We describe in this 
section the main tables characterising our sample. An ex- 
cerpt of each table is given here, while the complete tables 
can be found in the electronic edition of the paper. 

In Table[S]we list the coordinates, GALEX DR4, SDSS 
DR6 and UKIDSS DR4 magnitudes. Occasionally multiple 
SDSS and GALEX magnitudes are available for one system. 
In these cases we give averaged magnitudes. We used PSF 
(point spread function) SDSS magnitudes when available, 
fibre magnitudes otherwise. GALEX magnitudea are avail- 
able for 1327 WDMS binaries, UKIDSS magnitudes for 466. 

In Table|6] we provide relative numbers of the stellar 
components in our WDMS binary catalogue. For the white 
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Figure 6. Two-component fits to SDSSJ 082609.72+194126.3, a 
WDMS binary in our sample. The top panel shows the WDMS 
binary spectrum as black line, and the two templates, white dwarf 
and M-dwarf, as dotted lines. The bottom panel shows the resid- 
uals (in units of the standard deviation) from the fit. The SDSS 
spectrum identifiers MJD, PLT and FIB are given below the ob- 
ject name. See also Figs. 17181 



Table 6. The catalogue divided into the different WDMS binary 
subtype, as defined in Sec. [4] 



type 


number 


type 


number 


DA/M 


1176 


DB/M: 


1 


DB/M 


45 


DA:/K 


5 


DC/M 


34 


DA/K: 


10 


DA/K 


49 


DH/M 


l(x) 


DB/K 


3 


LARP 


1(*) 


WD/M 


135 


PG1159/M 


1(+) 


WD/K 


26 


(WD/M) 


32 


WD/M: 


1 


(WD/K) 


2 


DA:/M 


27 


(DA/M) 


4 


DB:/M 


1 


(DC/M) 


1 


DC:/M 


3 


(DA/K) 


1 


DA/M: 


43 







fx) ISzkodv et al l d2009h ; (*) ISchmidt et all 1120051) : (+) 
iNagel et al.l l|2006l ) 

dwarfs we use the flags DA, DB, DC, WD (if the white dwarf 
type is unknown), DH (if the white dwarf is magnetic), and 
PG1159 (hot hydrogen-deficient pre-white dwarf). We use 
the acronym LARP to indicate a low-accretion-rate polar. 
The secondary stars are flagged as M or K according to their 
spectral type. If the flag is followed by a colon, the classifi- 
cation of the stellar component is uncertain. Finally, we list 
in brackets the WDMS binary candidates. We consider basi- 
cally two types of candidates:(l) systems with very low S/N 
ratio spectroscopy that does not allow a clear classification 
of the stellar components in their spectra, and (2) systems 
with marginal blue (red) excess in their SDSS spectra, with- 
out morphological problems in their SDSS images ( Sec. 12.31) . 
and with \ 2 values (Sec. l2.2[ ) favouring a binary classification 
but no GALEX (UKIDSS) magnitudes available to confirm 
the existence of the second component. 



10 A. Rebassa-Mansergas et al. 



Table 5. The complete catalogue. Coordinates and GALEX-SDSS-UKIDSS magnitudes for the 1602 WDMS binaries and candidates 
are also included. The entire table (including also the photometric errors) is provided in the electronic edition of the paper. We use "-" 
to indicate that no magnitude is available. 



SDSS J 


ra[°l 

L J 


dec[°] 


nuv 


fuv 


U 


CI 


V 


i 


z 


V 


./ 


H 


K 


000152.09+000644.7 


0.46704 


0.11242 


18.45 


17.90 


19.03 


18.61 


17.94 


17.50 


17.25 


16.51 


16.05 


15.40 


15.28 


000442.00-002011.6 


1.17500 


-0.33656 






23.72 


20.38 


19.13 


18.65 


18.28 










000611.94+003446.5 


1.54975 


0.57958 


21.78 




21.38 


20.92 


20.12 


19.00 


18.38 


17.53 


17.05 


16.58 


16.20 


001029.87+003126.2 


2.62446 


0.52394 


20.17 


19.96 


21.92 


20.85 


19.97 


19.00 


18.42 


17.65 


17.14 


16.52 


16.36 


001247.18+001048.7 


3.19658 


0.18019 


20.50 


20.71 


20.73 


20.21 


19.66 


18.63 


17.96 


17.09 


16.60 


16.13 




001339.20+001924.3 


3.41333 


0.32342 


16.41 


19.73 


15.94 


15.56 


15.55 


15.63 


15.89 










001359.39-110838.6 


3.49749 


-11.14405 


17.77 


17.42 


18.30 


18.43 


18.31 


20.75 


22.82 










001549.02+010937.3 


3.95425 


1.16036 


20.97 


20.68 


21.23 


20.86 


20.60 


19.85 


19.27 


18.46 


17.86 


17.45 


17.12 


001726.63-002451.1 


4.36096 


-0.41419 


19.71 


20.30 


19.67 


19.28 


19.02 


18.18 


17.54 


16.60 


16.07 


15.56 




001733.59+004030.4 


4.38996 


0.67511 


20.83 


22.43 


22.09 


20.79 


19.58 


18.17 


17.38 


16.37 


15.84 


15.27 


14.97 


001749.24-000955.3 


4.45517 


-0.16536 


15.87 


15.40 


16.56 


16.86 


17.03 


16.78 


16.47 


15.75 


15.33 


14.76 


14.56 


001853.79+005021.5 


4.72412 


0.83931 


20.46 


20.27 


21.00 


20.38 


19.64 


18.80 


18.35 


17.52 


17.09 


16.51 


16.27 


001855.19+002134.5 


4.72996 


0.35958 


22.42 


22.12 


21.60 


20.60 


19.87 


18.97 


18.38 


17.54 


17.09 






002143.78-001507.9 


5.43242 


-0.25219 


22.20 




22.58 


19.63 


18.39 


17.02 


16.30 


15.40 


14.87 


14.31 


14.05 



5 STELLAR PARAMETERS 

In IRebassa-Mansergas et al] (l2007n : ISchreiber et al.1 (|2008h : 
iRebassa-Mansergas et alT i 20081 ) we presented a spectral de- 
composition/fitting technique and a M-dwarf spectral type 
- radius relation to determine the stellar parameters of 
WDMS binaries spectroscopically identified by the SDSS. 
Our method allows to estimate the effective temperature, 
surface gravity, mass and radius of the white dwarf, as well 
as the spectral type and radius of the main sequence com- 
panion. In addition, two independent distance estimates can 
be obtained by estimating the best-fit flux scaling factors o f 
the two components (see IRebassa-Mansergas et all (|2007T ) 
for details). We here basically use the same procedure but 
incorporate two modifications. 

Firstly, we compiled an additional library of 222 
high S/N spectra of D B white dwarfs from SDSS DR4 
|Eisenstein et al] 120061 ) covering the entire range of ob- 
served T c ff. Decomposing the binary spectrum and fit- 
ting the stellar compo n ents in the same way as in 
IRebassa-Mansergas et all l|2007l ) but using the just men- 
tioned DB templates instead of DA templates allows to es- 
timate the DB effective temperatures in our catalogue. 

As a second modification we take into account the ul- 
traviolet GALEX magnitudes if the spectral fitting does not 
provide a unique solution for the white dwarf temperatures. 
An example is given in Figs. 16181 Performing the spectral 
decomposition of SDSSJ 082609.72+194126.3 into the white 
dwarf and main sequence components (Fig.[SJ) and modelling 
the Balmer lines of the white dwarf model fit provides two 
solutions for the white dwarf temperature, the so called 
cold and hot solutions (black dots in the top right panel 
of Fig[7J . In most cases an additional fit to the entire spec- 
trum breaks this ambiguity and clearly indicates which of 
the two solutions has to be preferred. However, in the case 
of SDSSJ 082609.72+194126.3 the ambiguity remains as the 
fit to the entire spectrum results in values that fall exactly 
on the line of maximum H/3 equivalent width. As shown in 
the bottom panel of Fig. [8] the ultraviolet magnitudes mea- 
sured by GALEX clearly exclude the hot solution. 

The stellar parameters and distances obtained for the 



1602 WDMS binaries in our catalogue are given in Table]?] 
Both the cold and the hot solutions are provided for each 
spectrum, while the solution preferred by us is given in the 
first line. The complete table is available in the electronic 
version of the paper. 



6 DISTRIBUTION OF STELLAR 
PARAMETERS 

We present in this section distributions of surface gravity 
and effective temperature of the white dwarfs, as well as 
the spectral type distribution of the companion stars in 
our catalogue. To facilitate the comparison with previous 
works we also provide the distribution of white dwarf masses. 
Our template fitting algorithm covers secondary star spec- 
tral types M0-M9 and the spectral type distribution shown 
here includes only clear classifications in this range. For the 
distributions of the white dwarf stellar parameters we only 
considered systems with relative errors smaller than 25%. 
This resulted in 1433, 1198, 1127, and 558 WDMS binaries 
for the spectral type, surface gravity, effective temperature, 
and white dwarf mass distributions, respectively (see Fig. [9]). 

In gene ral terms the distributions are s imilar to those 
presented in IRebassa-Mansergas et al] (|2007l ): the most fre- 
quent white dwarf temperatures range between 10000- 
20000 K, white dwarf masses cluster around M w( j — 0.5Mq, 
loggr ~ 7.8 for the vast majority of the white dwarfs, and 
the spectral type of the companion stars are typically M3- 
4. In Fig. 1 101 we show the T c g, A/ W( j, loggr, and spectral 
type cumulative distributions obtained from the WDMS bi- 
naries s tudied in this work (blue lines) , the systems anal- 
ysed in IRebassa-Mansergas et ail (|2007l ) (red lines) and of 
a volume-limited sa mple of single white dwarfs (black lines, 
iHolberg et al] |2008). Kolmogorov-Smirnov (KS) tests were 
applied to compare the sets of effective temperatures, white 
dwarf masses, and logy. As the secondary spectral type dis- 
tribution consists of discrete values we performed a \ 2 test 
in this case. We briefly describe each parameter distribution 
in the following sub-sections. 
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Table 7. White dwarf masses, effective temperatures, surface gravities, spectral types and distances of the 1602 WDMS binaries in our 
catalogue, as determined from spectral modelling. We list both the hot and the cold solutions, with preferred solution given in the first 
line for each spectrum. The other solution is given for completeness. The complete table, including also notes for individual systems, 
can be found in the electronic e dition of the paper. We use t he flags e, s, a, and h an d "re " for those s y stems which have been studied 
previously by lEisenstein et al" | d2006h . ISilvestri etail [120071) , lAugusteiin et al.l [120081 1 and iHeller et all [|2009h , and which arc resolved 
WDMS binary pairs, respectively. Again, we indicate that no stellar parameters are measured with 



SDSS J 


type 


MJD 


PLT 


FIB 


Trff[k] 


err 


log g 


err 


MwdfM^l 

VVl_l L I , \ J 


err 


d w d [P c ] 


err 


Sp 


d scc [pc] 


err 


flag 


013356.07-091535.1 


DA/M 


53612 


1915 


431 


12392 


1324 


7.47 


0.36 


0.36 


0.16 


531 


112 


8 


217 


100 




013356.07-091535.1 


DA/M 


53612 


1915 


431 


12250 


1328 


7.49 


0.36 


0.36 


0.16 


52.3 


110 


8 


247 


100 




013418.52+010100.0 


DA/M 


53741 


1502 


517 


22811 


2601 


7.75 


0.41 


0.50 


0.22 


912 


224 


1 


1146 


225 


re 


013418.52+010100.0 


DA/M 


53741 


1502 


517 


9844 


473 


8.66 


0.46 


1.01 


0.25 


201 


73 


1 


1146 


225 




013441.30-092212.7 


DA/M 


52147 


662 


477 


12110 


2657 


7.06 


1.13 


0.24 


0.48 


2464 


1314 


2 


1914 


456 




013441.30-092212.7 


DA/M 


52147 


662 


477 


11699 


2808 


7.14 


1.19 


0.26 


0.56 


2315 


1285 


2 


1914 


456 




013441.30-092212.7 


WD/M 


52178 


662 


468 


















2 


1900 


453 




013441.30-092212.7 


WD/M 


52178 


662 


468 


















2 


1900 


453 




013504.31-085919.0 


DA/M 


53612 


1915 


464 


9187 


480 


8.82 


0.69 


1.10 


0.36 


218 


134 


4 


471 


139 


a 


013504.31-085919.0 


DA/M 


53612 


1915 


464 


25891 


4599 


7.40 


0.77 


0.38 


0.36 


2119 


919 


4 


471 


139 


a 


013716.08+000311.3 


DA/M 


52203 


698 


465 


18756 


3667 


8.52 


0.70 


0.94 


0.37 


586 


295 


2 


856 


204 


s/e 


013716.08+000311.3 


DA/M 


52203 


698 


465 


15601 


5326 


8.88 


0.95 


1.14 


0.49 


364 


326 


2 


856 


204 


s/e 


013750.03+003237.6 


DA/M 


52203 


698 


470 


22292 


894 


8.38 


0.14 


0.86 


0.08 


271 


29 


2 


272 


65 


s 


013750.03+003237.6 


DA/M 


52203 


698 


470 


16717 


1890 


9.24 


1.16 


1.32 


0.60 


91 


31 


2 


272 


65 


s 


013801.04+230329.3 
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6.1 White dwarf temperature 

We obtain a 40% probability (called KS probability in what 
follows) for the maximum vertical distance between the 
white dwarf te mperature distribution obtained here and the 
distribution in lRebassa-M anscrga s""et alj (|200?t ) being equal 
to or larger than the maximum vertical distance between 
the two cumulative white dwarf temperature distributions. 
Hence, there are no indications that the two distributions 
are not drawn form the same parent population. In contrast, 
comparing the WDMS binary white dwarf temperature dis- 
tribution with the temperature distrib ution of the volume- 
limited sample of single white dwarfs (|Holberg et al.l [20081 . 
black lines in Fig. llOp we obtain a KS probability of 10 -15 . 
This is straight forward to understand as the presence of the 
secondary star makes the identification of cold white dwarf 
primaries rather difficult. 

6.2 White dwarf mass and surface gravity 

Comparing the white dwarf mass and surface gravity dis- 
tributions obtained here with the sample described in 
iRebassa-Mansergas et all (|2007t ) we obtain KS probabilities 
of 50%, and 2% respectively. The low KS probability for the 
log g distributions is caused by the significantly higher frac- 
tion of low-mass white dwarfs in the present sample (see 
Sect l7.1fl . This is most likely an effect of the systematic 
changes in the flux calibration pipeline of SDSS spectra. 
These changes especially affect the blue part in the spec- 
trum and, hence, the fit to the continuum that allows to 
distinguish between the hot and cold solutions. In addition, 
the decision between hot and cold solution is affected by 



taking into account GALEX ultraviolet fluxes, one of the 
modifications incorporated in this paper (see also Sect. 17. 

Comparing the obtained distributions with those de- 
rived from a volume -limited sample of single white dwarfs 
<|Holberg et al l |200S| . black lines in Fig.rTUJ) we obtain KS 
probabilities of 10 -12 for the white dwarf masses and 10~ 15 
for the surface gravities. The white dwarf mass distributions 
of WDMS binaries and single white dwarfs are indeed ex- 
pected to be different. The white dwarf mass distribution 
of WDMS binaries is composed of two different types of 
systems, i.e. WDMS binaries with an initial main sequence 
binary separation large enough to avoid mass transfer dur- 
ing the evolution of the more massive star, and systems that 
suffer from CE evolution (PCEBs). The white dwarf mass 
distribution of the non-interacting WDMS binary popula- 
tion is expected to be identical to the mass distributions of 
single white dwarfs, i.e. clustering at M w( j ~0.6Mq, while 
PCEBs are expected to contain a large number of low- 
mas s He-core white dwarfs (M w d ^ 0.47Mq). For exam- 
ple, Ide Kool fc Ritteil (Il993l ) estimate the fraction of such 
systems among PCEBs to be about 50%. As ~ 30% of the 
presented sample of WDMS bi naries from the SDSS have 
evolved through a CE phase l|Schreiber et al.l [2008). the 
rather large fraction of white dwarfs with M w d — 0.4 Mq 
in Fig. is not surprising but expected. We are currently 
pursuing a large-scale follow-up project of WDMS binaries 
to identify the PCEBs in our sample and to measure their 
orbital parameters. A detailed analysis of the white dwarf 
mass distribution of PCEB and wide WDMS binaries is un- 
derway (Rebassa-Mansergas et al. 2009a, in preparation). 
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Figure 7. Spectral model fit to the white dwarf component of 
SDSSJ 082609.72+194126.3 obtained after subtracting the best- 
fit M-dwarf template. Top left panel: best-fit (black lines) to the 
normalised H/3 to He (gray lines, top to bottom) line profiles. 
Top right panels: 3, 5, and 10<r x 2 contour plots in the T e ff — log g 
plane. The black contours refer to the best line profile fit, the 
red contours to the fit of the whole spectrum. The dashed line 
indicates the occurrence of maximum H/3 equivalent width. The 
best "hot" and "cold" line profile solutions are indicated by black 
dots, the best fit to the whole spectrum is indicated by a red dot. 
Bottom panel: the residual white dwarf spectra resulting from 
the spectral decomposition and their flux errors (gray lines) along 
with the best-fit white dwarf model (black line) and the residuals 
of the fit (gray line, bottom). See also Figs l6l8l 
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Figure 8. Top: in blue the cold white dwarf solution obtained 
in Fig. [7] plus the best-fit M-dwarf template obtained in Fig. [6] in 
red the same but for the hot solution; in gray the SDSS spectrum 
of SDSSJ 082609.72+194126.3; in black dotted line the best-fit 
M-dwarf template, and the white dwarf models that satisfy the 
cold and hot solutions in Fig. [7] Bottom: the same but including 
the near- and far-ultraviolet GALEX magnitudes of this object 
(red dots). See also Figs l6l7l 
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Figure 9. White dwarf effective temperature, surface gravity and 
mass, and spectral type of the secondary distributions obtained 
from the SDSS WDMS binaries in our catalogue. 



6.3 Secondary star spectral type 

The distribution of spectral types is shown in the bottom 
panel of Fig. [9] The cut-off at early spectral types is very 
likely a consequence of selection effects as discussed in more 
detail in Sect .[5] The strongly decreasing number of WDMS 
binaries with late type (M7-M9) companions might also be 
related to selection effects, as late-type stars are faint and 
harder to detect against a moderately hot white dwarf. Nev- 
ertheless, SDSS covers a much broader colour space than 
previous surveys and, in principle, should be able to iden- 
tify more WDMS binaries containing cool white dwarfs plus 
very late-type com panions. It is worth mentioning also that 
iFarihi et all {2005) have constructed the relative distribu- 
tion of spectral types in the local M/L dwarf distribution, 
which peaks around M3-4, and steeply declines towards 
later spectral types, suggesting that late-type companions 
to white dwarfs are intrinsically rare (see Sect .[5] for more 
details). The \ 2 comparison of the WDMS binary spectral 
type distribution presented her e with the one described in 
iRebassa-Mansergas et all (120071 ) gives a probability of 82% 
that both distributions are drawn from the same parent dis- 
tribution (see also Fig. llOf) . 



6.4 Distances 

In Sect.[5]we derived two independent distance estimates for 
the WDMS binaries in our catalogue. We compare these val- 
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Figure 10. Effective temperature (top left), white dwarf mass 
(top right), surface gravity (bottom left) and spectral type (bot- 
tom right) cumulative distributions obtained from the stellar pa- 
rameters ij3f_the_Wj2MS_bjrmx^ here (blue lines) , anal- 
ysed in iReb assa-Manscrg as et all (120071) (red lines), and from 
a volume-limited s ample of single white dwarfs (black lines, 
iHolbere et aHl2008T) . 



ues in Fig. llll We here only consider systems with relative 
errors in the white dwarf distance less than 25%. The relative 
error in d sec is dominated by the scatter in the spectral type - 
radius relation provided in lRebassa-Mansergas et al.l (120071) . 
Hence, it represents an intrinsic uncertainty rather than a 
statistical error related to the fit and we therefore do not ap- 
ply any cut in d soc - Moreover, we excluded from the analysis 
any WDMS binary containing a white dwarf of less than 
12000 K, as spectral fits te nd to overestimate the surface 
gravity, and hence the mass IjKoester et al. 2009), leading to 
an underestimate of the white dwarf radius, and hence the 
distance. If more than one SDSS spectrum is available, we 
use averaged distances. This procedure resulted in a sample 
of 603 WDMS binaries. 

In the top panel of Fig. llll black dots represent systems 
in which the two distances agree within 1.5a (~3/4 of the 
total sample) level, while the differences in the two distances 
obtained for the objects in red exceed this limit. Most of the 
> 1.5cr outliers are above d scc = d w d (~l/5 of our WDMS 
binaries), while we find d scc < d w( j for only ~ 5% of the en- 
tire sample. Hence, there seems to be a systematic effect that 
leads to overestimating the seco ndary star distances. As in 
iRebassa-Mansergas et ail l|2007l ) we assume that magnetic 
activity raises the temperature of the inter-spot regions in 
active stars that are heavily covered by cool spots, leading to 
a blue-shift of the optical colours compared to inactive stars. 
However, to finally evaluate this interpretation one needs to 
perform a detailed activity analysis of our sample based on 
H a emission. This would further require to distinguish be- 
tween the different sub-samples forming the WDMS binary 
population, i.e. PCEBs and wide systems, as the fraction of 
active stars is expected to depend on the rotation rate. This 
is beyond the scope of this paper but will be presented in 
a forthcoming publication (Rebassa-Mansergas et al. 2009b, 
in preparation). 

As a first simple test, we adjust the spectral type of the 
secondaries to bring into agreement the two distances. It 



turns out that a change of 1-2 spectral subclasses is enough 
for the majority of cases. Only ten systems in the right panel 
of Fig.fTTJneed a change of more than two spectral subclasses 
to reach d scc — d w d. We have inspected these objects in 
more detail and the large discrepancy might be explained as 
follows: 

• Six objects contain hot white dwarfs, i.e. 
SDSSJ 003221.86+073934.4, SDSSJ 032510.84-011114.1, 
SDSSJ 080229.99+072858.1, SDSSJ 095719.25+234240.8, 
SDSSJ 101323.90+043946.1, SDSSJ 141536.40+011718.2. If 
these are short orbital period systems, irradiation of the 
secondary by the hot primary may lead to overestimating 
the distance to the secondary star. These WDMS binaries 
might therefore be considered as candidates to probe 
for radial velocity variations. SDSS J032510. 84-011114.1 
benefits from two SDSSJ radial velocity measurements in 
Tabled but no variation is detected. This may reflect the 
speculative nature of the just given argument or be caused 
by the SDSS spectroscopy sampling the same orbital phase 
twice. 

• Three systems (SDSSJ 025306.37+001329.2, 
SDSSJ 204729.04-064536.7, SDSSJ 210624.12+004030.2) 
are resolved in their SDSS images. Depending on the exact 
placement of the fibre, the flux contribution of one or 
both stars is likely to be underestimated. This translates 
to underestimated flux scaling factors and overestimated 
distances. 

• For one object, SDSSJ 232624.72-011327.2, the S/N of 
the corresponding SDSS spectrum is very low (S/N = 4.2). 
This is probably the reason for the discrepancy found in the 
distances. 



7 COMPARISON WITH PREVIOUS STUDIES 

Since ISmolcic et"al] (|2004T l discovered the WDMS binary 
bridge several SDSS WDMS binary catalogues have been 
presented, and several fitting routines to determine the stel- 
lar parameters have been applied. In this section we com- 
pare in more detail the results of our spectral decomposition 
method with those obtained by earlier studies. 



7.1 Comparison with Rebassa-Mansergas et al.l 
420071 ). 

Given that the S DSS spectra reduction pipeline was im- 
proved with DR6 l|Adelman-McCarthv et al.ll2008l ). we de- 
cided to compare the stellar parameters obtain e d her e with 
those we obtained in lRebassa-Mansergas et alj {2007). The 
stellar parameters for those objects with multiple SDSS 
spectra were averaged. Fig. [12] compares white dwarf effec- 
tive temperatures, surface gravities and secondary spectral 
types. Both studies agree within the errors in the major- 
ity of cases, with average relative differences of 13.5% and 
2.5% in T e ff and loggr respectively, and an average differ- 
ence of 0.3 spectral subtypes. We obtain significantly differ- 
ent values especially for systems containing white dwarfs 
with temperatures in the range 10.000-20.000 K (see top 
panels of Fig. ll2|l . In ~50% of the cases this is due to the 
fact that we are making use of ultraviolet GALEX mag- 
nitudes to constrain the white dwarf line fit solutions in 
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Figure 11. Left panel: Comparison of d sec and d w( j obtained from our spectral decomposition and white dwarf fits to the SDSS spectra. 
Approximately 1 /5 of the systems are dscc ^ ^wd outliers by more than 1.5(7 (red dots). Right panel: the spectral types of the secondary 
stars were adjusted by 1-2 spectral classes to achieve d wd = d aec . 



the present paper (see Sect. [5]). In the remaining ~50% 
this is most likely a consequence of the systematic changes 
in the flux calibration pipeline. The most dramatic case 
is SDSS J151045. 70+404827, a clear outlier in the bottom 
left panel of Fig. 1121 with a difference in T e g and logg of 
21.000 K and 1.7 dex respectively. This was one of the cases 
in which the solution was modified by the use of GALEX 
magnitudes. The obtained differences in spectral type are 
consistent with the general uncertainty of our decomposi- 
tion/fitting procedure of 0.5 spectral subclasses (see Sect. [5]). 
Only for SDSS J 173548.36+541424.4 this difference exceeds 
one spectral subtype. This system is very faint (i ~ 20), and 
the low S/N ratio is causing the spectral type determination 
to be rather uncertain. 

7.2 Comparison with Ivan den Be sselaar et al. 
(|2005h 

We have identified in this work 53 DB/M dwarf b i naries , 
13 of them in common with Ivan den Besselaar et al.l |2005). 
The stellar parameters of these systems are given in Ta- 
ble[8] Appa rently, the white dwarf tempera tures differ signif- 
icantly. As Ivan den Besselaar et al] (|2005h used DR3 SDSS 
spectra and measured all parameters by assuming a white 
dwarf mass of 0.6 M©, these discre pancies are not too sur- 
prisin g. In addition, as stated by Ivan den Besselaar et al.l 
l|2005l ). their derived effective temperature is related to the 
fitting of the secondary star, i.e. changing the spectral type 
by one subclass can lead to differences in the white dwarf 

effective temperatures of 8000 K to 10000 K. 

The distance measurements by Ivan den Besselaar et all 

(2005) are based on the white dwarf fitting, while we use 
the secondary star spectral type in case the primary is not 
a DA white dwarf. Taking into account the uncertainties in- 
volved in both distance measurements, we obtain reasonable 
agreement between both values. 

Concerning the spectral type of the compan- 
ion stars our values are in good agreem ent wit h 
those obtained by Ivan den Besselaar et al] {2005). 
In all but two (SDSS J093645. 14+420625.6 and 



SDSS J113609. 59+484318.9) cases the difference does 
not exceed one subclass. SDSS J093645. 14+420625.6 is 
a i = 20th magnitude object with a low S/N spectrum 
leading to large uncertainties in the obtained parameters. 
SDSS J113609. 59+484318.9 contains a hot white dwarf 
(>30000 K) that significantly contaminates the spectrum of 
the M-dwarf. 



7.3 Comparison with Isilvestri et"aTl (20061 

We compare in Fig.Q2] the white dwarf effective tempera- 
tures and surface gravities, and secondary spectral types 
obtain ed in Sect .[5] with those obtained by ISilvestri et al] 
(2006). We considered a sample of 421 spectra for which 
both studies present values for the stellar parameters. 

The average relative difference in effective tempera- 
tures and surface gravities is reasonably low, i.e. 14.7% 
and 2% respectively. For systems with white dwarf tem- 
peratures below ~20000 K, however, the obtained values 
can differ by up to 2210 K an d 1.6 dex. As discussed in 
iRebassa-Mansergas et all l|2007l ). we interprete this strong 
disagreement to be caused by the ambiguity between the hot 
and the cold solutions. At higher temperatures (> 50000 k) 
o ur fits tend to provide lower values of T e g than those given 
in lSilvestri et al] (|2006l ). This is probably caused by their use 
of DR4 spectra that were reduced with a different pipeline. 
In the majority of cases the secondary star spectral types 
are in reasonably good agreement (i.e. the difference is not 
exceeding one subtype). However, for ~17% of the WDMS 
binaries the difference is of two or more subtypes, with a 
maximum difference of four. 

W e have also inspected the systems that Silv estri et al] 
(2006) failed to fit and find that in ~70% of these we 
are able to find a solut i on. A s previously discussed in 
IRebassa-Mansergas et all l|2007h this indicates that our 
method is more robust if the S/N ratio is low or if one of 
the stellar components contributes little to the total flux. 
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Figure 12. Bottom panels from left to right: comparison of the white dwarf effective temperature s, surface gravities, and the spectral 
types of the secondary stars determined in this work and those of Rcbassa-Manscrg as et al, | d2007h . Top panels, from left to right: the 
white dwarf effective temperature ratio, and the difference in surface gravity and the secondary's spectral types from the two studies as 
a function of the white dwarf temperature. 
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Figure 13. Bottom panels from left to right: comparison of the whi te dwarf effect i ve tem peratures and surface gravities and the spectral 
types of the secondary stars determined from our fits and those of lSilvestri et al. (2006). Top panels, from left to right: the white dwarf 
effective temperature ratio, and the difference in surface gravity and the secondary's spectral types from the two studies as a function of 
the white dwarf temperature. 
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Figure 14. Bottom panels from left t o ri ght: comparison of the white dwarf effective temperatures, distances and surface gravities 
determined from our fits and those of iHeller et al.l «2009r) . Top panels, from left to right: the white dwarf effective temperature and 
distance ratios, and the difference in surface gravity from the two studies as a function of the white dwarf temperature. 
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Figure 15. Bottom pan els from left to righ t: comparison of the secondary star effective temperatures, distances and masses determined 
from our fits and those of IHeller et all {200g). The colors red, blue, green, yellow, magenta, black, and cyan refer to spectral types M1-M7 
respectively. Top panels, from left to right: the secondary star effective temperature and distance ratios, and the mass difference from 
the two studies as a function of the white dwarf temperature. 
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Table 8. Effect ive temperatures, spectral t ypes and distances 
obtained from Ivan den Bcssclaa r et al,l l|2005h and this work. The 
first line on each systems corre sponds to the results obtained by 
Ivan den Bcsscla ar et al,l d2005f) . the second line provides our re- 
sults. 
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7.4 Comparison with iHeller et all (|2009h 

iHeller et all (|2009l ) studied a sample of 857 WDMS bi- 
nary candidates. For 636 of these systems they provide 
the results of an independent spectroscopic parameter fit- 
ting method. FigureQ3] compares the white dwarf effective 
temperatures, dist ances and surface gravities obtained by 
IHeller et alj (|2009T ) with the values obtained here. The aver- 
age relative differences between the two analysis are 12.1% 
(white dwarf temperature), 23.9% (distance) and 3.6% (sur- 
face gravity) . Th e agreement betwe en both studies is dif- 
ficult to asses as IHeller et al.l (l2009t ) do not provide error 
estimates. If one assumes uncertainties of 2000 K and 0.3 
dex for their T e s and loggr values, and a typical error of 
17% for their obtained distances, we find that only 18%, 7% 
and 16% of the given values (for T e B, distance, and logg) 
do not overlap at the 1.5 a level. Hence we conclude that 
the results obtained in both studies are in agreement. Fi- 
nally, we note that the horizontal patt erns in the bottom 
right panel of Figure [H illustrate that IHeller et ail (120091 ') 
estimate the surface gravities by using a grid with a rather 
poor resolution of 0.5 dex. 

Fig. [15] compares the effective temperatures, distances 
and masses of the secondary stars in both data sets. In this 
case the average relative differences are 3.5% (effective tem- 
perature), 42.1% (distance) and 4 2% (mass). While w e used 
M star templates in our fitting, IHeller et all (l2009h used 
a gri d of PHOENIX model spectra (Hauschil dt fc Baronl 
1 19991 ). The colors refer to objects containing secondaries of 
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Figure 16. Selection effects in SDSS WDMS binaries can be un- 
derstood by analysing the density maps obtained from their stel- 
lar parameters. From top to bottom the (log T c g,d wd ), (d wd ,Sp), 
and (log T c ff ,Sp) density maps. 



the same spectral type, i.e. red, blue, green, yellow, magenta, 
black, and cyan for spectral types M1-M7. Values of effective 
temperatures and masses for the systems in our catalogue 
were estimated using the spectral type-radius-mass-effective 
tempe rature relations given in iRebassa-Mansergas et all 
(2007). Therefore WDMS binaries containing secondaries 
of the same spectral type are associated to the same ef- 
fec tive tempera t ure an d mass, while the estimates provided 
bv IHeller et all l|2009h cover a considerably larger range in 
both mass and effective temperature (lower left and lower 
right panels of Fig. [15}. According to the lower right panel, 
Heller et all ((2009) measure systematically lower masses. A 
comparison of the obtained distances to the secondary stars 
is shown in the middle panels of Fig. 1151 Since the observed 
flux is the same for a given SDSS spectrum, the only dif- 
ference between both methods comes from the radius and 
the flux at the stellar surface. As the grid of model spectra 
is not finer than our grid of template spectra both methods 
differ by a constant for a given spectral type/secondary mass 
combination. This can clearly be seen in the bottom middle 
panel of Fig. 1151 for each pair of template/model spectrum 
we obtain one or two straight lines in the dista nce-distance 
diagra m. It seems that we overestimate and /or IHeller et all 
(2009) underestimate the di s tance s. However, we presented 



Rebassa-Mansergas et all (|2007h (see also Sect. 16. 4[) a rea- 
sonable explanation for the ~l/5 of the WDMS binaries 
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with d scc > d w d, while there is no obvious physical mecha- 
nism that may account for the estimated white dwarf dis- 
tances being systematic ally smaller than t he companion dis- 
tances as the values of l|Heller et alj|2009h seem to suggest. 
We have to keep in mind though the uncertainties related to 
both methods are quite large and the activity interpretation 
of the systematically larger secondary distances still needs 
to be confirmed. We therefore conclude that the question of 
whether our templates or the model spectra provide more 
reliable distances remains open. 



8 SELECTION EFFECTS 

In lRebassa-Mansergas et al.l ([2007) we briefly discussed pos- 
sible selection effects that may affect the observed SDSS 
WDMS binary population. The much larger sample pre- 
sented here allows to investigate these selection effects in 
more detail. To avoid contaminations from unreliable stel- 
lar parameters we base our analysis on systems with mean 
relative errors of less than 25% in their white dwarf pa- 
ramet ers. As discussed in detail in lRebassa-Mansergas et al.l 
(|2007h the distances derived from the white dwarf parame- 
ters are probably more reliable than those derived from the 
secondary stars. We therefore use the white dwarf distances 
here and quote them simply as distances in the following. 

Fie. 1161 shows three density maps that illustrate the se- 
lection effects affecting the SDSS WDMS sample. Due to 
the restrictions in the white dwarf parameters 597, 692, and 
1052 systems (from top to bottom) were considered. The 
log T c ff , d w( j density map in the top panel shows that binaries 
with white dwarf primaries cooler than 10 000 K are only de- 
tectable at relatively short distances (d<400pc), while sys- 
tems containing hotter white dwarfs have been detected at a 
much wider range of distances (~ 200—1000 pc). In addition 
there is obviously a general trend of increasing distance with 
white dwarf temperature. This is straight forward to under- 
stand as cold white dwarfs become too faint to be detected 
at larger distances. In contrast hot white dwarfs are intrin- 
sically brighter, but also rarer than cold white dwarfs, and 
hence dominate at larger distances where the volume sur- 
veyed by SDSS is sufficiently large. Most objects are hence 
concentrated at ^400-500 pc, with white dwarf effective tem- 
peratures between ^15000-25000 K, as also shown in the ef- 
fective temperature distribution (Fig.0. The middle panel 
of Fig. 1 161 shows a similar effect but for the secondary star 
spectral types, Sp. Early M-dwarfs are relatively hot and 
consequently exceed the lower SDSS brightness limit at dis- 
tances of < 300 pc. In contrast, late- type secondaries are cold 
enough to be detected at short distances (~ 100—200 pc) but 
too dim to be observed at distances larger than ~ 500 pc. 
Finally, the bottom panel of Fig. 1161 shows the \ogT B s, Sp 
density map. A clear trend of later spectral type compan- 
ions to colder white dwarf primaries can be seen. This is 
again easy to understand: while late-type companions to hot 
white dwarfs are too faint to be detected, cold white dwarf 
primaries are out-shined by early spectral type secondaries. 

The selection effects just described can explain the cut- 
off at early spectral types in the bottom panel of Fig[9] as the 
white dwarf primaries are not detectable. The scarcity of sys- 
tems with later-type (>M6) secondaries, however, is proba- 
bly not only related to selections effects, as the spectral type 



distribution of low-mass field stars also peaks at Sp ~ M4- 
5, and decline towards later s pectral types fe.g. lFarihi et alj 
|2005| ; iReid et al l [20071 . |200S|) . The lack of WDMS binaries 
with late-type companions is therefore probably an intrin- 
sic property of the WDMS binary population that appears 
more pronounced due to selection effects. We will address 
this systematically in a forthcoming publication. 

From the analysis of Fig.[l6]we conclude that a "typical" 
SDSS WDMS binary contains a M3-4 companion, a ~10000- 
20000 K primary, and is observed at a distance ~400-500 pc. 
However, we have to keep in mind that a typical SDSS 
WDMS binary is not necessarily a typical WDMS binary. 
Overcoming the selection effects just described requires to 
combine the SDSS with complementary magnitude limited 
surveys. Detecting binaries consisting of a hot white dwarf 
and a late-type companion is most likely to arise from the 
use of infrared surveys such us UKIDSS or 2MASS, while 
the identification of cool white dwarfs with early-type dom- 
inated M-dwarfs requires to incorpor ate ultraviolet surveys 
such as GALEX l|Maxted et al.ll2009h . 



9 COLOUR-COLOUR DIAGRAMS 

We have provided in the previous sections of this paper a 
detailed spectroscopic analysis of WDMS binaries in SDSS. 
In this section we make use of the photometric magnitudes 
given in Tableland combine them with the stellar param- 
eters measured from the SDSS spectra to investigate the 
appearance of WDMS binaries in colour-colour space. Fig- 
ure !17l and [TH] show four relevant colour-colour diagrams. 
Stellar sources are shown in gray, DA/M WDMS binaries 
in yellow, the few DB/M are in blue and the DC/M binaries 
in green. Finally, DA-DB-DC/K systems are shown in red, 
while the black dots represent the WD/M-K binaries (see 
Sect. l2~4)l . 

A general feature evident in all diagrams is that a cer- 
tain number of systems appear to be outliers from the gen- 
eral stellar locus of WD MS binaries, i.e . the W DMS binary 
bridge described first in ISmolcic et al.l (|2004l ). We have in- 
spected these outliers and find the majority of them being 
resolved in the SDSS images (see Sec. l2.3|l . 

To further evaluate the information provided by pho- 
tometry we show the white dwarf effective temperature and 
secondary star spectral type distributions in six different 
regions for each colour-colour diagram, as indicated by the 
horizontal (dotted) and vertical (dashed) lines. The resulting 
distributions are shown above and below the colour-colour 
diagrams. As previously, we considered only those white 
dwarf effective temperatures with relative error less than 
25%. In the following sub-sections we briefly describe the 
main conclusions that can be drawn from the four colour- 
colour diagrams. 



9.1 u - g vs g - r 

The most commonly used SDSS colour-colour diagram for 
stellar sources in SPSS is u — q ys q — r |Fanl Il999l ; 
iRichards et~aH |2002| ; ISchreiber et ail 120071 ). Inspecting the 
left hand side of Fig.[l7]it can clearly be seen that the white 
dwarf effective temperature distribution is shifted towards 
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Figure 17. Middle left: u — g vs g — r colour-colour diagram. WDMS binaries are represented according to their binary components 
as follows: DA/M binaries in yellow, DB/M in blue, DC/M in green, DA-DB-DC/K in red, and WD/M in black. Stellar sources are 
represented with gray dots. Two vertical (dashed) lines divide the diagram in three rectangular regions (columns left, centre, right). In 
the same way, two horizontal (dotted) lines divide the diagram in three different regions (rows top, middle, bottom). Top and bottom 
left: distributions of white dwarf effective temperature and spectral type of the companions obtained for the six different regions outlined 
above (three rows and three columns). The right panels follow the same structure, but for the nuv — u vs u — r colour-colour diagram. 



lower temperatures if u — g or g — r increases, i.e. for red- 
der colors (compare the left, centre, and right as well as the 
top and the middle distributions) . Only two WDMS binaries 
with reliable effective temperatures are found in the bottom 
region close to the main sequence. As expected, the distri- 
butions of secondary spectral types contain more early-type 
secondaries if one moves to redder colors (left to right or top 
to bottom). 

9.2 nuv - u vs u -r 

The nuv — u vs u — r colour-colour diagram is provided in 
the middle right panel of Fig. 1171 The white dwarf effective 
temperature and spectral type distributions in the six re- 
gions of the diagram (see top and bottom right panels in 



Fig. ll7|l are similar to those discussed in the above subsec- 
tion, the main difference being the decrease in the number 
of systems. This is due to the smaller percentage of SDSS 
WDMS binaries that have been detected with GALEX. 

9.3 nuv - i vs nuv - H 

Figure[l8] (left panels) shows the nuv — i vs nuv — H colour- 
colour diagram and the corresponding distributions of white 
dwarf effective temperatures and secondary spectral types. 
Clearly, requesting SDSS-GALEX-UKIDSS magnitudes re- 
duces the number of WDMS binaries and the shown dis- 
tributions are statistically less robust. However, the gen- 
eral trend observed in the previous colour-colour diagrams 
remains: systems composed of white dwarfs hotter than 
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Figure 18. Same as in Fig. 1171 but for the nuv — i vs nuv — H (left) and i — J vs J — H (right) colour-colour diagrams. 



20000 K are generally detected in the top and left regions 
of the diagram and there is a clear trend of decreasing white 
dwarf effective temperature towards redder colours (top to 
bottom and left to right). The previously observed shift to- 
wards earlier spectral types for redder colours seems still to 
be present but is much less pronounced. The most striking 
feature of the colour-colour diagram is the nice separation 
of WDMS binaries and single stars. 



9.4 i - J vs J - H 

Finally, we provide in Fig. [18] (right panels) the i — J vs 
J — H red colour-colour diagram and the corresponding dis- 
tributions. Again, due to the reduced number of systems 
the overall trend of having less hot white dwarfs and more 
early spectral type secondaries for redder colours seems to 
be present but less significant. Nearly all systems contain- 
ing hot white dwarfs are located in the upper region. This 
colour-colour diagram represents an additional example of 



nicely separating WDMS binaries and single stars. We pro- 
vide several colour-cuts that can be used to select WDMS 
binaries in the next section. 



10 COLOUR CUTS 

Having studied in the previous section the relation between 
colours and stellar parameters (i.e. white dwarf effective 
temperature and spectral type) in four colour-colour space 
diagrams, we define here four colour-cuts of WDMS bina- 
ries. The fact that we are considering a total of 11 photo- 
metric band passes (two from GALEX, five from SDSS and 
four from UKIDSS) increases considerably the number of 
possible colour-cut selections. We here provide four exam- 
ples (two of them already introduced in Sect .[9]) in which 
WDMS binaries are clearly separated from the locus of sin- 
gle main sequence stars (see Fig. llDJl . To quantify how com- 
plete these colour-cuts are would imply an analysis on the 
different SDSS-I, SEGUE, Stripe 82 areas, where different 
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Figure 19. Colour-cut selection for WDMS binaries in nuv — i vs nuv — H (top left), i — J vs J — H (top right), fuv — nuv vs r — z 
(bottom left), and y — K vs u — z (bottom right). Main sequence stars and white dwarfs are represented in light gray, WDMS binaries 
in dark gray, and colour-cuts by black straight lines. 



target strategies were tested. Such an endeavour is beyond 
the scope of this paper and will be pursued elsewhere. 

The top left panel of Fig.ll9lshows the nuv—i vs nuv — H 
colour-colour diagram introduced above, which offers an ex- 
cellent opportunity to unambiguously isolate WDMS bina- 
ries from single main sequence stars and white dwarfs with a 
simple colour-cut. We represent in light gray main sequence 
stars and white dwarfs, in dark gray WDMS binaries. With 
a straight black line, i.e. 



(nuv - i) < -0.85 + 0.83 x (r 



H) 



(2) 



both populations can be distinguished. 

As we have seen in Sect.|9j WDMS binaries and sin- 
gle stars also separate nicely in the i — J vs J — H colour 
plane. However, the location of WDMS binaries containing 
hot white dwarfs overlaps with those of single stars. The 
colour-cuts shown in the right top panel of Fig. 1191 i.e. 



0.3 < (J - H) < 0.7 
(i — J) > 1.2 
J) > 0.4 + 1.85 x( J-H) 



(3) 
(4) 
(5) 



will therefore mainly select WDMS binaries composed of 

cold white d warfs. 

Recently iBianchi et al.l (|2007T ) studied the properties of 
the GALEX-SDSS matched source catalogues and classified 
sources by studying their colours. Inspired by their Fig. 5, 
we show in the left bottom panel of Fig. [19] the fuv — nuv 
vsr-z colour diagram, and provide colour-cuts that should 



select the main population of WDMS binaries: 

(r - z) > -0.3 (6) 
(fuv - nuv) < 0.85 + 3.9 x (r - z) (7) 

Finally, we provide in the bottom right panel of Fig. 1191 
the colour-colour diagram and colour-cuts for WDMS bi- 
naries in y — K vs u — z. Again stars are represented in 
light gray, WDMS binaries in dark gray, and colour cuts by 
straight black lines . This colour diagram has been already 
used bv lChiu et all (|2007T ) for quasar selection. We here only 
slightly modified their colour-cuts and obtain 

(u - z) > 5 (8) 
(y-K)> 0.7 (9) 
(y - K) > 0.6 + 0.14 x (u - z) (10) 

which should successfully select WDMS binaries. 



11 RADIAL VELOCITIES AND NEW PCEB 
CANDIDATES 

In this final section o f th e paper we follow 
iRebassa-Mansergas et al.l l|2007h and use the 
NaiAA 8183.27,8194.81 absorption doublet and/or 
the Hcv emission to measure radial velocities. As in 
IRebassa-Mansergas et all (120081 ); ISchreiber et all i|2008h we 
use a single width parameter for both line components in 
the Nai doublet. The radial velocities of 1068 systems with 
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Table 9. Radial velocities measured from the 
Na I AA 8183.27,8194.81 doublet and the Ha emission for 
1068 systems in our catalogue. The complete table can be found 
in the electronic edition of the paper. In the last column we 
quote with "y" and "n" those radial velocity values obtained 
from spectra that are, and are not combined from individual 
exposures taken on different nights, respectively. We use "-" to 
indicate that no RV is available. 



Table 10. Upper limits to the orbital periods of the five PCEB 
candidates identified in Sect. 1111 White dwarf masses are taken 
from Table[7] except for SDSS J2346+4340, where we assume of 
mass of 0.5 M«^n (see Fig-EJ) - Secondary star masses are estimated 

from Table 5 in Rebassa-Manserg as et al ] (|2007tl . K sec values are 
obtained from Table[9] where we use the Nal doublet radial ve- 
locities for SDSS J2346+4340. 
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pronounced Nai AA 8183.27,8194.81 absorption doublet 
and/or Ha emission in their SDSS spectra are given in 
Table|9] where we also include the HJD of the o bserva - 
tiomjjT As discussed in iRebassa-Mansergas et all l|2007| ). 
a significant fraction of SDSS spectra are combined from 
observations that have been taken in different nights which 
leads to a reduced sensitivity to radial velocity variations. 
We indicate this in the last column of Table[9] 

A comparison with t he ra dial velocities obtained in 
IRebassa-Mansergas et ail l|2007| ) gives an average relative 
difference of 19.5%, and the measurements generally over- 
lap within the errors. The reasons for the small changes have 
been a lready discussed in detail in lRebassa-M anscrga s~et al.l 
(2008) and can be summarised as follows: (1) we modified 
the procedure to fit the Nai absorption doublet by using a 
single width parameter for both line components; and (2) 
we used DR6 spectra here (instead of DR5) . 

As described in the introduction, the orbital period dis- 
tribution of WDMS binaries is expected to be bimodal, sep- 
arating long orbital period systems whose stellar compo- 
nents evolve like single stars from short orbital period sys- 
tems that suffe red from mass transfer interactions, mostly 
CE evolution llWillems fc Kolbl I2004T ). As demonstrated 
by IRebassa-Mansergas et al] ( 20071 ). multiple SDSS spec- 
troscopy of WDMS binaries can be used to measure radial 
velocity variations thereby eventually identifying candidates 
for being a PCEB. Our criterion for calling a WDMS bi- 
nary a strong PCEB candidate are radial velocity variations 
with 3cr significance. This is evaluated by calculating the \ 2 
of the measured radial velocities against their mean value. 



SDSS J 


0333(*) 


0743 


1453 


2318 2346 




+0054 


+2835 


+0055 


+0034 +4340 


Porb[d] 


575 


121 


8 


5 2 



2 Note that the HJDs in IRebassa-Manserg as et al. I n »Tt i were 
wrong by -0.5 days because of an erroneous conversion of the FITS 
headers of the SDSS spectra. 



(*) The effective temperature is below 12000 K, and consequently 
the white dwarf mass is likely overestimated. 

The incomplete 7 function then gives the probability for the 
measured radial velocities being consistent with a constant 
value. If this probability is below 0.0027 the measured ra- 
dial velocity variations are 3er significant and we consider the 
corresponding WDMS binary a strong PCEB candidate. 

We find here nine and 16 PCEB candidates using the 
Nai doublet and Ha emis sion respectively. A comparison 
with the results obtained in lRebassa-Mansergas et al] (|2007l ) 
shows that four systems, i. e. SDSS J030904. 82-010100. 9, 
SDSS J113800.35-001144.5, SDSS J173727.27+540352.2, 
and SDSS J234534.50- 001453.7 are not anymore considered 
PCEBs, while we find five new PCEB candidates, namely 
SDSS J033301.51+005418.5, SDSS J074329.62+283528.0, 
SDSS J145300.99+005557.1 SDSS J231874.73+003403.3, 
and SDSSJ234638.76+434041.7. Apparently, the dif- 
ferences in the method of determining radial velocities 
and/or the re-reduction of the SDSS data can cause a 
given system to move e i ther way across our criterion 
l|Rebassa-Mansergas et all 120081 ). The Nai absorption 
lines from the secondary star are a more robust probe for 
radial velocity variations, and we consequently consider 
that additional follow-up spectroscopy is necessary for 
those PCEB candidates identified from Ha radial velocity 
variations. Upper limits to the new PCEB candidate orbital 
pe riods have been estimated in th e same way as described 
in lRebassa-Mansergas et al] {2007). These values are given 
in TablefTOl 



12 SUMMARY 

We have presented a catalogue of 1602 WDMS binaries from 
the spectroscopic SDSS DR6. We have used a decomposi- 
tion/fitting technique to measure the effective temperatures, 
surface gravities, masses and distances to the white dwarfs, 
as well as the spectral types and distances to the companions 
in our catalogue. Distributions and density maps obtained 
from these stellar parameters have been used to study both 
the general properties and the selection effects of WDMS 
binaries in SDSS. A comparison between the distances mea- 
sured to the white dwarfs and the main sequence compan- 
ions shows d scc > d w( j for ~l/5 of the systems. We also have 
made use of GALEX, SDSS and UKIDSS magnitudes to 
study the distribution of WDMS binaries in colour-colour 
space and present simple colour-cuts that allow to clearly 
separate WDMS binaries from other stellar objects. Finally, 
we have measured radial velocities for 1068 WDMS binaries 
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measured from the NaiAA 8183.27,8194.81 absorption dou- 
blet and/or the Ha emission line. Among the systems with 
multiple SDSS spectroscopy, we find five new WDMS bina- 
ries showing significant radial velocity variations identifying 
them as PCEB candidates. The here presented new, up- 
dated, and most complete catalogue of WDMS binaries from 
the SDSS represents a superb data base for future follow-up 
studies that may significantly contribute to a better under- 
standing of close compact binary star evolution. 
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